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The dispersion of a coupled resonator optical waveguide (CROW) made of photonic crystal mode- 
gap cavities is pronouncedly asymmetric. This asymmetry cannot be explained by the standard tight 
binding model. We show that the fundamental cause of the asymmetric dispersion is the inherent 
dispersive cavity mode profile, i.e., the mode wave function depends on the driving frequency, not the 
eigenfrequency. This occurs because the photonic crystal cavity resonances do not form a complete 
set. We formulate a dispersive mode coupling model that accurately describes the asymmetric 
dispersion without introducing any new free parameters. 


In coupled cavity systems in photonic crystals (PhC), 
light transport occurs via evanescent field coupling be¬ 
tween high-quality factor (Q) cavities. The well known 
example is the coupled resonator optical waveguide 
(CROW) which is a linear chain of cavities [T]. Arrays of 
coupled cavities have attracted substantial scientific at¬ 
tention for practical applications such as slow light engi¬ 
neering and strong light-matter interaction enhancement 
US], and many novel phenomena of fundamental inter¬ 
est such as gauge fields and time-reversal of light 

pulses [Hill- 

Realizations of low loss and compact CROWs require 
high-Q, wavelength-sized cavities. Mode-gap cavities [5] 
created by shifting some of the holes of around PhC 
waveguides as shown in Fig. 1(a) have been demon¬ 
strated to be extremely suitable for creating large scale 
cavity arrays nnuni. 

The tight binding (TB) model [U [12] is the usual ap¬ 
proach for modeling the dispersion in coupled cavities 
[311 dun]. The core concept of the TB model is that 
in coupled cavities the wave functions are tightly confined 
in each cavity. The eigenmodes of the individual cavities 
are then coupled to yield waveguide modes. Coupling 
between neighboring cavities is due to the overlap of the 
eigen wave function of the cavity modes. The dispersion 
of the CROW band as predicted by the TB model is a 
cosine curve [T]. For a given structure, its parameters can 
be evaluated by numerical methods, such as plane wave 
expansion (PWF) |15j and finite difference time domain 
(FDTD) [11 . 

However, the TB model does not describe the disper¬ 
sion of a CROW composed of mode-gap cavities correctly. 
In Fig. 1(b) we show the group velocity curve as ob¬ 
tained from the TB model and as calculated by FDTD 
m- The TB result is perfectly symmetric with respect 
to k = O.^tt/R. In contrast, the FDTD result is pro¬ 
nouncedly asymmetric. The asymmetric spectrum has 
also been observed experimentally HU. It is remarkable 
that the TB model fails here, as it only depends on the 



FIG. 1. (a) Schematic of a CROW composed of coupled mode- 
gap cavities. The cavities are created by shifting of the yel¬ 
low, red and purple holes around the waveguide in a tapered 
way. The width of the waveguide is 0.98-\/3u {a is the lattice 
constant). The radius of the holes is 0.25 a. The shift of 
the yellow holes is Sl=77 0.0124%/3a, is a factor we call it 
modulation strength. The shifts of the red and purple holes 
are 2/3 SI and 1/3 SI respectively.(b) Group velocity of the 
2D CROW consisting of mode gap cavities with modulation 
strength 1.5 calculated by 2D FDTD (black dots) and TB 
model (green line). Parameters are calculated by FDTD. 


assumption of energy-independent eigenmodes, which is 
a direct consequence of the completeness of the set of 
eigenmodes Cl¬ 
in this Letter, we uncover the physical origin of the dis¬ 
crepancy and propose an improved mode coupling model. 

The resonant modes of a closed conservative cavity 
form a complete set. An important consequence of com¬ 
pleteness is the fact that the response of a cavity is de- 
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fined in terms of modes do not themselves depend on the 
driving frequency. In an open cavity we have quasimodes 
which form a complete set only if certain conditions are 
fulfilled [12], essentially, the edge of the resonant mode 
must be clearly defined in the structure and no outgoing 
waves should be scattered back to the cavity. These con¬ 
ditions are never fulfilled in a PhC cavity. In a PhC cavity 
light is confined by constructive interference of Bragg re¬ 
flection and the Bragg reflection takes place throughout 
the PhC structure. As a result, the quasimodes do not 
form a complete set and their spatial profile depends on 
the driving frequency, i.e, mode profile is dispersive. 

The breakdown of completeness for an open cavity 
shows a signature in driven oscillation. If the cavity is 
driven at frequency w, the field inside the cavity will be 
E(r, t) = E(r,a;) exp(—where E(r, w) ^ E(r,a;o) . 
Alternatively speaking, when an open cavity is driven, 
the wave function of the mode is determined by the driv¬ 
ing frequency not the free oscillation frequency. 

In CROWs, therefore, the coupling between neighbor¬ 
ing cavities is not characterized by the overlap of the 
eigen modes of the single cavity E(r, wq) (wq is the eigen 
frequency of the cavity) but of the dispersive modes 
E(r,w). 

In Fig. 2 the dispersion of a quasimode is shown qual¬ 
itatively. In a true eigen mode, the size of the mode 
profile should be always the same as the cavity is driven 
at different frequencies. However, the quasimode is dis¬ 
persive and will be spatial narrower (Fig. 2(a)) or wider 
(Fig. 2(b)) than the resonant mode when the driving 
frequency is lower or higher than the intrinsic frequency 
respectively. 

-None dispersive 



FIG. 2. In driven oscillation, the quasimode (red) will be 
smaller (a) or larger (c) compared to the intrinsic reso¬ 
nant mode in the case when driving frequency is lower or 
higher than the intrinsic frequency. In contrast a true eigen¬ 
mode (black) keeps the same spatial profile when driven off- 
resonance. 

We performed two dimensional (2D) FDTD simula¬ 
tions to quantitatively study the dispersion of the mode 
profile of PhC mode-gap cavities. In the simulation a 
continuous source is placed in the center of a single cav¬ 
ity and switched on smoothly. After the transient died 
out, we output the Ey field across the waveguide direc¬ 
tion. The same procedure was repeated several times at 
different frequencies and different cavities [20] . 

In Fig. 3(a) we plot the amplitudes of the Ey field 




FIG. 3. (a) 2D FDTD simulation of forced driving oscillation 
of a mode-gap cavity. The amplitude of the Ey field across 
the waveguide direction versus the position is plotted. The 
center position is located at 0. The red, black and green lines 
illustrate cases with different driving frequencies (tJo=0.22212 
c/a, (5i=0.00011 c/a and 52=0.00010 c/a), (b) Red crosses 
(blue squares) depict decay lengths of a cavity with intrin¬ 
sic frequency a;i=0.22212 c/a (a;2=0.22222 c/a). The black 
curve is the model we use to describe the decay length. The 
decay lengths are obtained by fitting the envelope of the field 
from a; = 4a to a; = 14a. The only difference between the 
geometries of the two cavities is the modulation strength. 


versus the position across the waveguide direction when 
we drive the cavity at three different frequencies. When 
the driving frequency is wq + i5i which is higher than the 
resonant frequency lcq, the decay in space is slower. Al¬ 
ternatively speaking, the mode is larger in space. When 
the driving frequency is ujo — ^ 2 , the mode is smaller. This 
clearly shows that for a single mode-gap cavity, the mode 
profile differs depending on the driving frequency. Out¬ 
side the modulation part of the cavity, the envelope of 
the mode profile decays exponentially. Along the waveg¬ 
uide direction x, the envelope decays as exp(— x/Zd), 
where Id is the decay length and x = 0 is the cen¬ 
ter of the cavity. The decay length is related to the 
dispersion of the waveguide where the cavity lies in. 
The waveguide dispersion above the bandedge can be 
expanded as w = Wedge + {k — / a/)^ /(2m) [2T], where 

m = (9^w/9^/c)“^ is the photon mass in analogy to the 
effective mass of electrons m and Wedge is the frequency 
of the edge of the waveguide band. By analytic continua- 
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tion of the dispersion, the wave vector k becomes a com¬ 
plex number when w < Wedge ■ The decay length which 
is the inverse of the imaginary part of the wave vector, 
follows as Id = {2m{u}edge — In Fig. 3(b), we 

plot the calculated decay lengths by the analytical model 
and the extracted decay lengths from FDTD simulations 
of two different cavities at various driving frequencies. 
Firstly we see when the cavity is driven at different fre¬ 
quencies the decay lengths differ. Secondly for two dif¬ 
ferent cavities (one with intrinsic frequency Wi and the 
other with intrinsic frequency 102 ), when the driving fre¬ 
quencies are the same, the decay lengths are the same. 
The decay lengths are well described by the analytical 
model, and we conclude that mode-gap cavities have a 
mode prohle that depends on the driving frequency not 
the resonant frequency. 

As the consequence of the dispersive mode prohle, the 
mode functions in the expression of the coupling rate 
[Dua should be the dispersive mode (DM) not the eigen 
mode. The expression of the coupling rate becomes 

a;/(5£(r-R)E,^(r-R) ■E^(r)dr 
/e(r)E<^(r) • E,^(r)dr 

In Eq.(l), E^(r) is the electric held of the wave function 
at frequency w, e(r) is the dielectric constant of a single 
cavity and 5e(r — R) is the dielectric difference between 
one isolated cavity and two cavities with heart to heart 
distance R. This expression is identical to that derived 
by Haus et aZ.[22] and by Yariv et al. [1], except for the 
replacement of the resonant mode prohle E^^ (r) by the 
dispersive mode prohle E^(r). The resulting dispersion 
of the CROW is, 

w = Wo-I-A-t-F(a;) cos(/ci?), (2) 

where wq is the eigen frequency of the single resonators, 
A is the frequency diherence between wq and the center 
of the CROW band. 



FIG. 4. Normalized coupling rate r(w)/ro of the CROW 
structure in Fig. 1 versus the normalized frequency detuning 
5. 


To evaluate Eq.(2) and (3), some assumptions are 
made. Firstly, we assume that each polarization of E,^(r) 
can be decomposed into {x)E^^ (jj, z), here x, y repre¬ 
sent the spatial coordinates and x is along the waveguide 
direction, where the E^^{yz) can be approximated to be 
non dispersive, E^^{y, z) « E^^{y, z). Secondly, we ap¬ 
proximate the envelop of E^ (x) by 1/ cosh(a;/Z_D) , which 
is a resonable approximation for mode-gap cavities. The 
last assumption is for 5e(r — R) and e(r). We divide 
(5e(r — R) into 5 blocks. For each block, the effects of Se 
is essentially the increase of the effective dielectric func¬ 
tion of the waveguide. We approximate that this increase 

happens at the center of each block, so we express Se as a 
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collection of delta functions Se = ^ S{x — Xj)Sej with 

i =-2 

Xj = ja and Sej = Seo (where Seo is a constant). 
Based on these three assumptions the coupling rate is 
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r(w) = pY1 

t =-2 


(3-|j|)^ 

Id cosh(A^)cosh(g) 


(3) 


Here /3 is a proportionality constant resulting from the 
integral over y, and we may abbreviate r(w) = /3g(w) 
where we emphasize g{u}) can be evaluated analytically. 
The dispersion relation becomes 


w = Wo -I- A -1- (ig{ijj) cos(fci?), (4) 

there are only three unknowns (w. A, /3) which is exactly 
the same number of the unknowns the TB model has, i.e., 
our model does not introduce any extra free parameters. 

We show the normalized coupling rate of the CROW 
structure in Fig. 1 as a function of normalized frequency 
detuning 5 = (w —Wo)/ro withFg = r(wo) in Fig. 4. The 
normalized coupling rate is dehned as r(w)/Fo, it is also 
equivalent to g{uj)/g(uJo). The normalized coupling rate 
in Fig. 4 increases non-linearly as the detuning increases, 
in other words, F(w) is highly dispersive. At the low 
k values which correspond to positive detuning S, the 
coupling rate is larger, and therefore the group velocity 
is enhanced. As a result, the maximum of the group 
velocity shifts to low k values as compared to the result 
from the TB model. 

We present the calculated group velocity of the CROW 
in 2D and in a 3D membrane structure in Fig. 5. In Fig. 
5(a) the results of group velocity calculated by DM model 
is asymmetric with a maximum at fc « 0.10 (2'k/R), this 
matches the 2D FDTD results extremely well. In Fig. 
5(b) the group velocity from the DM model is asymmet¬ 
ric with maximum at A: « 0.13 (2tt/R), it agrees satisfac¬ 
torily with the fully 3D FDTD results |23]. In contrast, 
the result from the TB model is symmetric about k = 0.5 
{tt/R). The excellent match between our DM model and 
the FDTD data confirms that dispersion of the cavity 
quasimodes is indeed the physical reason for the failure 
of the TB model to accurately describe photonic crystal 
CROWs. 
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FIG. 5. (a) Group velocity of the 2D GROW consisting of 
mode gap cavities with modulation strength 1.5 calculated 
by FDTD (black dots), TB model (green line) and DM model 
(red line). All the parameters are as same as in simulations 
of driven oscillation, (b) Group velocity of the 3D GROW 
consisting of mode gap cavities with modulation strength 2 
calculated by FDTD (black dots), TB model (green line) and 
DM model (red line). In the 3D simulation, h = 0.37a (h is 
the thickness of the slab) and e = 10.04. 


In conclusion, the dispersion of a CROW composed 
of photonic crystal cavities can not be described by the 
standard TB mode due to the breakdown of complete¬ 
ness of the resonant modes. We show a new model tak¬ 
ing into account of the dispersive property of the mode 
profile of the single cavity in a CROW structure. Our 
model describes the dispersion of the PhC CROW accu¬ 
rately without additional free parameters. Dispersion of 
the modes is inherent in all cavities. For racetrack res¬ 
onators nm] and PhC cavities with resonances far away 
from the bandedge such as L3 cavities [53], dispersion is 
negligible and shows no discernible effect on the disper¬ 
sion of the related CROW structures. However, for shal¬ 
low defect cavities, such as mode-gap cavities with small 
modulation strength and double-heterostructure cavities 
PSIITT] with small lattice mismatch, the defect modes are 
close to the edge of the waveguide band which make the 
dispersive nature of the modes become pronounced. The 
dispersion of CROWs of such cavities has strong asym¬ 
metry. Thus, our theory will be useful for describing all 
devices based on coupling of shallow defect cavities such 
as delay lines m. optical memory |28j . and VT symmet¬ 
ric diodes 
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